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abstract 

As pan of the National Aero-Space Plane ("AST) 

effects of gravitational force initial tank FLOW- 
heat transfer rate on the two-dimensional ‘empire ^ ed for the evaluation. These 

3D, a commerdal ffnhe-^fference fluid flow ^ M experimental data with 

effects were examined on the bas , P ^ ^ com pared against an existing 

gaseous hydrogen P“f^“‘ddftion the effect, of mesh sine and convergence criteria 
one-dimensional model. In addi , „ u foI f uture modifications and 

on the analytical results were investigated. Suggest oils lor 
uses of FLOW-3D for modeling a NASP tank are a so presell e 

INTRODUCTION 

The temperature distribution in the ullage 

has been of interest for design in ot gr °“" ‘ . j^cant in determining the energy 

terns. The ullage gas temperature is n , . t \ xe ramp pressurization process 

distributions and pressurant gas requirem j^-4 p rev ious experimental and 

as well as during the pressurized expulsion ofthetan^ ^ J resu l ts which 

analytical studies on tank pressunza ion a ture distri b u tions and pressur- 

show the effects of various parame ers on jj-io jj owe ver, these studies were 

ant gas quantities required during 1 * 8 ® ^lodeUng, and the gravitational forces were 
generally constrained to one-dimens ' k the National Aero-Space 

fimited to those found at ground-level ^ te ftt va nous gravita- 

Plane (NASP) includes the study of f ^^J““ herm ; dynaini c effects. As part 
tional levels with the potential for mu . ij d and liquid hydrogen, is also 

propellant "i^ 1 seen^uring^ liquid hydrogen propellant. Therefore, a study was mitiated 
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to determine the multi-dimensional effects on the temperature profiles in the ullage 

gas during the ramp pressurization and expulsion phases on the operation of a NASP 
vehicle tank. 

The FL0W-3D computer code was chosen for modeling the multi-dimensional 
temperature profiles which occur during tank pressurization. FLOW-3D is a gen- 
eral three-dimensional fluid flow program developed by Flow Science, Inc. 13 - 14 . The 
code uses a finite-difference solution of the Navier-Stokes equations based on the 
M ar ker- an d - Cell (MAC) technique and the Implicit-Continuous Fluid-Eulerian (ICE) 
technique 16 ’ 16 . These techniques were developed through a series of codes with the 
generic name SOLA. SOLA-VOF, the most successful of these codes, serves as the ba- 
sis for FL0W-3D; details of SOLA-VOF can be found in reference 17. Capabilities of 
LOW-3D include incompressible or fully compressible flow, rotating and accelerating 
reference frames, and turbulent transport modeling. FLOW-3D has been used exten- 
sively to model tank systems, especially systems in which sloshing is involved 18-20 . 


FLOW-3D CODE SIMULATIONS 


The 1988 version of FL0W-3D was used in the present study to determine the 
temperature profiles during the pressurization phase of operation of a tank. In the 
pressurization process the tank pressure is gradually increased by the addition of a 
gas through a diffuser. There is no liquid outflow during pressurization. In a previ- 
ous experimental study 7 a 5 foot diameter spherical tank filled with liquid hydrogen 
propellant was pressurized from approximately 17.4 psia to 50 psia using gaseous 
hydrogen. The tank geometry from this study serves as the basis for the FLOW-3D 
modeling discussed here. As part of the NASP project, this 5 ft tank geometry will 
also be used in tests to be performed at K-Site, NASA Lewis Research Center’s slush 
hydrogen test facility located at Plum Brook Station in Sandusky, Ohio 12 . 

A two-dimensional representation of the 5 foot tank was used in the FLOW-3D 
simulations with ullage volume percentages of 55%, 28%, and 8.5%, as shown in Figure 
1. The simulated ullage was actually a very thin slice taken from the center of the 
tank; however, because the thickness of the slice was small (0.1 ft), the representation 
was essentially two-dimensional. The ullage fractions used correspond to actual liquid 
hydrogen test data discussed in reference 7. In this simulation, hydrogen gas enters the 
ullage space through a hemispherical diffuser, modeled as a porous obstacle (an option 
within FLOW-3D). The computational grid used for most cases was a rectangular 10 
cell by 10 cell mesh. The wall of the tank cut through this mesh, and the cell size 
remained constant throughout the mesh for all cases studied. 
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A W tVip<;p simulations was determined from 
The velocity of the entering gas use ^ the smaller simulated ullage 

data in reference 7. This v ^ locl y was SC h d f t h e pressurization phase 

volume into account such that the ^ 0 f P temperal „res were 

was near 50 psia, correspond, ng to the P«™us da a 1 g 2g4 . R for 

also modeled from the , S 5% ullage Pressurization times used in the 

the 28% Ullage, and 524” R for the 8.5% «£ge , d for most cases the initial 

computations were oMarne^ lhe casc w here the ramp pressurization 

— . ^ ^ Bec 2 rse R in a stahi.iti: in the 

tank, decreasing linear y was uged with the incompressible liquid, only 

code surfaced when a compres g . , , /.l- i s known as an “impedance 

the compressible gas in the ^ oMd and the gas is 

sb trrr zxzz * — - 

these simulations, as well as some significant outputs. 

For the three ullage volume ' ^ ^ 

eters on the ullage tempera ure pr initial ullage temperature (T;), ramp 

gravity level (G), initial tank pressure ( .)» , heat trans f er rates. In the 

pressurization profile, and wall (Qwa ) compared against results 

heat ^^g^^e^^Res^ndi Center one-dimensional ramp pres- 

^"T^tL, results are presented showing effects of convergence 

criteria (e) and mesh sire on FLOW-3D results. 


results 


Gravity 

ro To 3 an d 4 show the effect of variations in the gravitational forces on 
Figures 2 , 3 , and 4 snow an d 8 5% ullage fraction, respectively. The 

temperature contours for 55%, o, ^ indicated above each graph, 

highest and lowest contour temperature^ ^ ^ presented cn each graph ; 

and correspond to H and L on tl g P ronstan t difference in temperature, as 

the distance between the contours represei includes the convergence 

deffned by the h,gh and low the analys.s. Each case 

cntena and time step ( e an ’ P ia a „d fina | press ure of 50 psia, 

was run using an initial tank pressure oi P 
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corresponding to the datn in the literature. These cases were run assuming no heat 
transfer to the wall or hqmd. The liquid level is indicated by the bottom of the 
rame. The graph in the upper left-hand corner of Figure 2, at 1-g (G=32.2 ft/sec 2 ) 
shows that the final temperature contours (temperatures at the end of pressurization) 

♦“MSn flat ’ WltH a tem P erature contour «nge of 376*R near the top of the tank 
to 109°R near the gas/liquid interface. This indicates a temperature profile which 
is only axially- dependent, with no radial dependence at 1-g. As the gravitational 
forces are decreased the temperature profiles begin to show some radial dependence 
However, at 55% ullage fraction this dependence does not become significant until 
the gravitational force is less than 0.1 g (3.22 ft/sec 2 ). At 0-g the profiles are highly 
dependent on radius. In addition, the profiles appear to be unsymmetric at 0-g 
The unsymmetnc profiles may be due to instabilities in the code at low gravitational 
levels; these instabilities will be discussed later in the report. Note that at 0-g the 
assumption of a rigid plane at the liquid surface probably does not apply, and a more 
rigorous analysis would be required in this region. 

One possible measure of where gravitational forces become significant is the di- 
mensionless Froude number, Fr, which is defined as 


where u is the gas velocity (ft/sec), g is the gravitational force (ft/sec 2 ), and A is a 
characteristic length (ft). The Froude number represents a comparison between the 
inertial forces and the gravitational forces; the inertial forces do not become greater 
an the gravitational forces until the Froude number is greater than 1. For the tank 
geometry here, with a flow velocity of approximately 0.5 ft/sec and a characteristic 
length corresponding to the mesh cell length (approximately 0.25 ft), Fr greater than 
1 corresponds to a gravitational force less than 0.1-g. Therefore, the temperature 

profiles are not expected to show a significant change in shape until the force due to 
gravity is less than 0.1-g. 


. ulla ® e ) shows a similar result to that obtained with a 55% ullage 

At 8.5% ullage shown in Figure 4, the final temperature profiles become dependent 
on radius at a higher gravitational level (between 0.1 and 0.5-g). In these runs the 
high temperature contour is actual greater than the inlet gas temperature. This 
result was obtained due to adiabatic compression in the tank. In addition, the 8.5% 
ullage cases appear to exhibit larger differences in high and low temperature contours 
than those of the 55% or 28% runs, and some unsymmetrical profiles appear at the 
lower gravitational forces. Code instabilities related to the convergence criteria may 
be the cause of the unsymmetrical profiles; a discussion of the convergence criteria 
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- >. ", ,r™" sr-fi“ 

XnTnt^glound should apply to the NASP vehicle. 

Tnitial Pressure 

In cases where normal boiling 

initial pressure, prior to the ramp press ^ pressure may be as low as 1.1 

In the case of slush hydrogen, owe , Tliere £ ore) ^ ecause s l us h hydrogen may 
psia, the triple point pressure o y r g ■ ^ ’ re on temperature profile 

be used to fuel the NASP vehicle, the (fleet of imti^ P>* ^ profiles at the end 

was studied- Figure^) shows * comp^ with no heat transfer at the wall or liquid, 
of pressurization with a 55 /o g 6 na l temperature profiles 

In the top figure, at an however, the fluid is 

show no radial dependence. an 313°R), and the temperature 

well mixed (the temperature contours range The mixing of the 

shows a significant radial depen enc compar i S on is made between the velocity 

ullage gas can be seen m Figure , w e above the graphs show the scale for 

profiles at the two initial pressures (enu ^ 1? 4 pgia the ve locity of the gas 

velocity vectors, in ft sec). At an in P P velocity is significant 

to characterize this phenomenon. 

Similar results are seen when the ullage 373°R to 

17.4 psia initial pressure, the temperature . q-tej*" * , , ; a sUow 
109°R, as shown in ranged 330 to 

significant mixing and radial depe” d ' sures can be se en in Figure 8, where 

306”R. The difference in mixing at t P ^ penetrate thro ugh the ullage to 

velocity profiles are compare • hydrogen pressurant remains near 

the liquid surface at 1.1 P«a, 17 , 4 P™ *“ 7 comparisons at 8.5% ullage, 

the diffuser. Finally, Figure 9 shows the temp the lower pressure, leading 

Here, as in the previous cases, the mixing - enhanced at the 1 ^ p ^ ^ ^ 

to less stratification The ‘em P e-t”re S may be due to some 
case do appear to be unsymmetnc y .^Pj pressures . However, the 

ZSs^ cd" Appears to be a strong function of the initial tank 

pressure, regardless of ullage fraction. 
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Initial Ullage Temperature 

the T' ysis The t graph 8hows ,he finai io«“ 

£££r. 7£z alters n th ft ^ 

i*-- 1 r r:r 

throughout tfc- 

emperature is high, but the temperature throughout the rest of the t t • , C 
In contrast, when a temperature range of 250°R to 40° R is chosen at the sta W 

V*'t a 1 Mtly linear ’ rangmg fr ° ra to 

Ramp Profiles 

a^SSSS:- 

=sS=~SS3 

„ tho Pre8S “ re ’ aS n d,SCUSSed in « fere ”“ ?■ figure 11 shows the two profiles “ed 

=1 tty ZT ProfitT n't 1 " r C ° n8tanl S| ° P ' Pr ° fila ' » constant ^ 

pre::;ia«::trFig„rtr.w: tt w «:r M : in ih t* rata “-«■*- *• 

LTXt “L 8 gur Vth e 8h r of the 

resulting^ temperature^aT ^e^ea^of^iessarization^ tIILT t ° f “* 

Zz^V'jr lhal « 8 the - ttzttzz 

C;.r^ -pit tr etf ‘ h ;,~‘r r 


6 



Hpat. Transfe r 

For the previous cures examined it "tlZ ^ Z - 

the wall or liquid tromttM gas. n 1 jn , he determination of the stratifies- 

the heat transfer to the liqui p y , t transfer on temperature contours 

;;rii— ^ (one -“ nal rep ~ 

resentation) were investigated. 

Heat transfer results are shown in Kgure^ 

figure three cases are considered^ (upper right graph), and 

(185 BTU/hr) from the gas to the wall and to b J om grap h). 

80 ft-lb/sec (370 BTU/hr) to the ^ “ d " ^ith SLURP, the NASA Lewis 
The heat transfer rates were chosen ase ^ amic parameters during slush 

Research Center code for calcula ing SLUR p provided heat transfer rates to 

rrr^^ - — “ iun5 

based on area. 

The results of the heat transfer comparison 

the heat transfer does not significant y a included however, the magnitude of 

With heat transfer to the wall and the liquid The adi abatic 

the final temperatures of the ullage does ecrea the tank 

solution shows a temperature re contour, using 80 

” tht 

F o, comparison of the effects iff heat 

additional runs were made. Figure 14 shows the insults of the P flom 

no heat transfer to the waU or liquid (top graph^o ft ft/- l ^ ^ ^ 60 

the gas to the wall and to the liquid (center graph), an J run for M 

ft-lb/sec (278 BTU/hr) to the liquid (bottom graph). Th™ ra(es ^ 

Sen Wd on mns with SLURP, as discussed in the previous figure. 

From the output in Figure 14 it can be “'"^V^^Hhe 

« «W-= ^ ‘h= Hqmd surface (bottom 
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figure) the profiles did show some radial dependence, but the degree of dependence 
is not high. This radial dependence may be a result of some code instabilities as 
will be discussed m the following figures. As expected, as the heat transfer rate was 
increased the temperature of the ullage gas decreased. Under adiabatic conditions 
the temperature contours ranged from 627°R to 149°R (top figure). The temperature 
is higher than the inlet gas temperature due to adiabatic compression. At 60 ft-lb/sec 
to the liquid the ullage gas temperature ranged from 553°R to 108°R. Therefore as 
discussed in the previous section on initial temperature differences, the heat transfer 
at the interface and axial ullage temperature profiles are related. 


igure 15 shows the comparison of centerline temperatures (the one-dimensional 
solution) for FLOW-3D runs with and without heat transfer. These results are com- 
pared against a similar run with SLURP. For these cases an ullage of 8.5% and an 
initial pressure of 17.4 psi were used. From the figure it can be seen that the adia- 
baticsolution gave temperatures up to 150°R higher than the solution obtained by 
SLURP If however a heat transfer to the wall and liquid is chosen such that it is 

rimwln calculated by SLURP, similar temperature distributions are calculated 
by r LUW-oD. 


, I o£° te t ^ at thC hCat transfer rates to wal1 a nd to the liquid were input 

to the FLOW- 3D code for study of the tank pressurization problems presented here. 

In these runs heat transfer to the wall could not be calculated as a mesh large enough 
to be used to model a volume composed of both a thin wall region and a relatively 
large ullage space would have been computationally prohibitive. In addition the 
impedance mismatch” between a low density hydrogen gas and a relatively high 
density liquid hydrogen propellant caused instabilities in the code. This prevented 
analysis of the heat transfer to the liquid from the gas, and limited the ability to 
analyze expulsion problems. Therefore, modifications to the basic FLOW-3D code, or 
a different approach to the problem (such as using FLOW-3D’s cylindrical coordinate 
capability) may be required in order to accurately predict temperature profiles in 

a NASP vehicle tank without prior experimental knowledge of the mass and heat 
transfer. 


Con ve rgen ce Criteria 

The following figures show the effect of convergence criteria on temperature con- 
our results obtained using FLOW-3D. The simulations were performed using a 55% 
ullage and an initial pressure of 17.4 psi, with no heat transfer to the wall or the liquid. 
The comparisons of various convergence criteria were examined both at 1-g and at 

lrr^rl y A 1Cal as3 ' riiptote of °-«- The convergence criterion, e (known as EPSI in the 
h LOW-3D input), is used to determine the termination of iterations - the iteration 
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continues until every cell in the computational mesh converges to within this value of 
e In the 1988 version of FL0W-3D the time step, dt (DELT in the input) and the 
convergence criterion are related because of stability limits on the code. Therefore, 
as the convergence criterion is decreased the time step must also be decrease . 

Figure 16 shows the effect of convergence criteria at 1-g. From the figure it can 
be seen that for the various convergence criteria used, .02, .002, and x 10 e 

is little difference in the results. The ullage temperature contours range 376 

at the top of the tank to 110-120°R near the gas/liquid interface. Although the 
temperature profile near the midpoint of the ullage differs somewhat between he 
cases considered, the convergence criteria does not seem to significantly affect the 
temperature distribution results at 1-g. This result is confirmed by a comparison 
the centerline temperatures, as shown in Figure 17. Note however that in some cases, 
such as those shown for the 8.5% ullage fraction, a loose convergence criteria (large 
value of “e”) may lead to code instabilities for the geometry chosen. 

Figure 18 shows the effect of convergence criteria at 0-g. In this figure the re- 
sults at a convergence criteria of 1x10-' is quite different from the results at e-0^ 
FVom the figure it appears that the temperature contours are more symmetrical a 
e=lxl0-' compared to the result with e=0.02. The symmetrical temperature dis- 
tribution is a result of an “even" distribution of the gas- The difference m veloaty 
profiles, as shown in Figure 19, indicates that the gae >' ^enly distributed when 
convergence criterion of 1*' is chosen (lower graph). At larger values of e th s 
symmetrical distribution of the gas is not evident. A symmetrical tempera ure ^d. - 
tribution is probably a realistic result for the geometry considered here ' th " et ° r ’ 
is recommended that a tight convergence criterion be used at low gravitation 
els (a much more computationally costly requirement) and that a loose convergence 
criterion be used for most cases at 1-g. 

Mesh Size 

The following figures show the effect of mesh size on temperature contours _ The 
results shown were for 55% ullage at 1-g and zero-g with an initial pressure of 17.4 
psi pressure convergence criterion (e) of .02 psi, and a time step (dt) of .005 sec. 
There was no heat transfer to the wall or liquid in these cases. Two mesh sizes were 
considered: a 100 cell mesh (10x10) and a 400 cell mesh (20x20). From the result 
shown in Figure 20 the mesh size did not impact the shape or the magnitude of the 
final temperature profiles to any large degree, although there was some difference 
near the liquid surface. From results presented in the Figure 21, however, mesh , 
appears to have a significant effect on temperature distribution at zero-g. I. 
case the 10x10 mesh may not be large enough to resolve the gas velocity profil 
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the tank. Therefore if a low gravitational force is used a finer partitioning of the 
tank volume should be used (at the cost of high CPU time). At ground conditions 
however, a smaller number of mesh cells appears adequate. 


FUTURE FLOW-3D MODIFICATIONS 

In order to use the FLOW-3D code for a complete analysis of the NASP vehicle 
ank several modifications to the code should be implemented. The ability to model 
e heat and mass transfer at the gas/liquid interface appears to be a critical capa- 
bility in determining the gas requirements, the slush hydrogen solid losses, and the 
egree of thermal stratification in the ullage 3 . The addition of this capability may 
require modification of the code to include the ability to perform calculations during 
expulsion, when the fluid dynamics of the liquid and the thermodynamics of the gas 
are both important. In addition, because the NASP vehicle may use both helium 
an ydrogen gas during the pressurization and expulsion phases, a two-component 
ulkgu model should be added. Finally, an initial slush hydrogen fluid dynamic 
model has been added to the FL0W-3D code 14 . However, this model will require 
improvements to fully model slush hydrogen dynamics, including an improved slush 
e ing/sohdification capability, the ability to have a free surface (gas/liquid inter- 
face) present when the slush model is used, and the capability of modeling slush 
hydrogen with a gaseous pressurant (3 phases simultaneously). 


CONCLUDING REMARKS 


The FL0W-3D computer code was used to determine the effects of various param- 
eters on the two-dimensional temperature profiles during ramp pressurization prior 
o pressurized expulsion of a liquid hydrogen tank. These results apply to the NASP 

vehicle tanks where multi-dimensional effects may be important. The results obtained 
can be summarized as follows: 


. For l-g and for an initial pressure near atmospheric, the ullage gas final tem- 
atures do not show a radial dependence. As the gravitational forces are reduced, 
w a significant radial dependence until the gravitational forces 


are less than 0.1-g. 


2. Initial tank pressure has a significant impact on the temperature distribution 
t an initial pressure of 17.4 psi the final ullage temperature profile was not dependent 
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on radius. At an initial pressure of 1.1 psi the final temperatures indicated that the 
fluid was well-mixed, and the profiles appeared to be radially dependen . 

3. The initial temperature affects only the magnitude of the final temperatures m 
the ullage. 

4. The shape of the ramp pressurization profile did not significantly impact the 
final temperature profiles. 

5 Convergence criteria and mesh size seem to make the largest difference in re- 
tional instabilities. 


6 Centerline temperatures using the SLURP one-dimensional pressumat, on code 

and the FL0W-3D program matched when heat transfer rates obtained by SLURP 
were scaled and used as input to the FL0W-3D simulations. 

As the current effort was intended to be an initial investigation of FLOW-® 
further analysis is required to determine the 

to analyze tank pressurization and expulsion pro • FLOW-3D such as the 

thermodynamics, FL0W-3D has been shown to be a versatile tool 
used for future studies of the NASP storage and feed system. 
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APPENDIX: Sample Input Files For FLOW-3D Runs 


Example Input File for a 552 Ullage Case 

K-SITE LH2 : 552 Ullage, G=32.2,Pi=17.4,Ti=250, 

$XPUT 

ITB=0 , EPSI=0.02, IPDIS=0, DELT=0.005, DTMAX=0.05, 
SPRTDT=4.0, PRTDT=5.0, PLTDT=4.0, 

VB=2, WT=6, WF=1 , WBK=1 , VL=1, WR=1, 

TWFIN=24.0, GZ=-32.2, 

IPUN=10, IHTC=1 , IFENRG=2, 

AUTOT=0 . 0 , 

TIMBCT(1)=0.0, TIMBCT(2)=0.2, TIMBCT(3)=100.0, 

TBCT(1,6)=307 .0, TBCT(2,6)=307.0, TBCT(3,6)=307.0, 
TBCT( 1,5) =40.0, TBCT(2,5)=40.0, TBCT(3,5)=40.0, 

TBCD=307. , 

WBCT(1,6)=0.0, WBCT(2,6)=-0.60, WBCT(3,6)=-0.60, 

PBCT(1,6)=3600. , PBCT(2,6)=3600. , PBCT(3,6)=3600.0, 
FBCT(1,6)=0.0, FBCT(2,6)=0.0, FBCT(3, 6)=0.0, 

P0BCT(l,5)=-40. , P0BCT(2,5)=-40. , P0BCT(3,5)=-40. , 
ICMPRS=1 , IADIX=1 , IADIZ=1 , 

RF2=24690. 96, CV2=53335.0, 

SEND 

$MESH 

PX(l)=-2.5, PX(2)=0.0, PX(3)=2.5, 

SIZEX(2)=0.25, 

NXCELL(1)=5, NXCELL(2)=5, NXCELT=10, 

PY(1)=0.0, PY(2)=0. 10, NYCELT=1 , 

PZ(l)=-0. 187 , PZ(2)=2.5, NZCELT=10, 

SEND 
$0BS 
N0BS=2 , 

CC(1)=6. 25, CX2(1)=-1.0, CZ2(1)=-1 .0, 

I0H(1)=1 , RAL(1)=0.5, 

CC(2)=6.0, CX2(2)=1 . 0, CZ2(2)=1.0, CZ(2)=-5.0, 

OP0R(2)=0. 7, I0H(2)=1 , 

TOBS(1)=0.0, TOBS(2)=100.0, 

P0BS(1, l)=-40. ,P0BS(2, l)=-40. , 

SEND 

$FL 

PRESI=2500.0, 

SEND 

$BF 

SEND 

$TEMP 

NTMP=1 , ITDIS( 1)=1 , 

TCC(l)=-52.82, TCZ=121 . 1 , TZH(1)=2.5, TZL(1)=0.75, 
SEND 
SMOTN 
SEND 
SGRAFIC 

NVPLTS=1, JV1(1)=2 , JV2(1)=2, 

NCPLTS=1 , K0NTYP(1)=5, JC1(1)=2, JC2(1)=2, 

SEND 

SPARTS 

SEND 
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Example Input File for a 28% Ullage Case 


K-SITE LH2: 28% Ullage, G=32.2,Pi=17.4,Ti=250, 
■XPUT 

ITB=0, EPSI=0.007, IPDIS=0, DELT=0.005, DTMAX=0.5, 
SPRTDT=4.0, PRTDT=5.0, PLTDT=4.0, 

WB=2, WT=6, WF=1 , VBK=1, WL=1, VR=1, 

TWFIN=24.0, GZ=-32.2, 

IPUN=10, IHTC=1 , IFENRG=2 , 


AUTOT=0.0, 

TIMBCT ( 1 ) =0 . 0 , 
TBCT(1,6)=284.0, 
TBCT(1,5)=40.0, 
TBCD=307 . , 
WBCT(1,6)=0.0, 
PBCT(1,6)=3600. , 
FBCT(1,6)=0.0, 
ICMPRS=1 , 
RF2=24690.96 , 


TIMBCT(2)=0. 2 , TIMBCT(3)=100.0, 
TBCT(2,6)=284.0, TBCT(3,6)=284.0, 
TBCT(2,5)=40.0, TBCT(3,5)=40.0, 


WBCT(2 , 6)=-0. 323 , WBCT(3, 6)=-0. 323, 
PBCT(2,6)=3600. , PBCT(3,6)=3600.0, 
FBCT(2,6)=0.0, FBCT(3 , 6)=0. 0, 
IADIX=1 , IADIZ=1 , 

CV2=51732.0, 


$END 

$MESH 

PX(l)=-2.5, 
SIZEX(2)=0.25, 
NXCELL(1)=5 , 


PY(1)=0.0, 
PZ(1)=0.75, 
SEND 
SOBS 
N0BS=2, 
CC(1)=6.25, 
I0H(1)=1, 
RAL(1)=0.5, 
CC(2)=6.0, 
0P0R(2)=0.7 , 
T0BS(1)=0.0, 


PX(2)=0.0, 

NXCELL(2)=5, 
PY(2)=0. 10, 
PZ(2)=2.5, 


CX2(1)=-1.0, 


PX(3)=2.5, 

NXCELT=10, 
NYCELT-1, 
NZCELT=10 , 


CZ2(1)=-1.0, 


CX2(2)=1.0, CZ2(2)=1.0, 

I0H(2)=1 , 

TOBS(2)=100.0, 


SEND 

$FL 

PRESI=2500.0, 


CZ(2)=-5.0, 


SEND 

$BF 

SEND 

$TEMP 

TCC(l)=-52.82, TCZsl21.lt TZH(1)=2.5, TZL(1)=0.75, 

SEND 

$M0TN 

SEND 

SGRAFIC 

NVPLTS=1, JV1(1)=2, JV2(1)=2, 

NCPLTS=1 , K0NTYP( 1 )=5 , 

JC1(1)=2, JC2(1)=2, 

SEND 

SPARTS 

SEND 
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Example Input File for an 8.52 Ullage Case 


K-SITE LH2: 8.52 Ullage, G=0.0,Pi=17.4,Ti=250, 
$XPUT 


ITB=0, EPSI=0.05, IPDIS=0, DELT=0.003, 
SPRTDT=4.0, PRTDT=5. 0, PLTDT=4.0, 

WB=2, VT=6, WF=1, WBK=1, VL=1, WR-1, 
TffFIN=15.0, GZ=0.0, 

IPUN=10, IHTC=1 , IFENRG=2, 

AUTOT=0.0, 


DTMAX=0 . 05 , 


TIMBCT(1)=0.0, 

TBCT(1,6)=524.0, 


TBCT( 1,5) =40.0, 
TBCD=307. , 

WBCT ( 1 , 6 ) =0 . 0 , 
PBCT( 1,6) =3600. , 
FBCT(1,6)=0.0, 
ICHPRS=1 , 
RF2=24690.96, 

$END 

$MESH 


TIMBCT(2)=0. 2 , TIMBCT(3)=100.0, 

TBCT(2,6)=524.0, TBCT(3,6)=524.0, 
TBCT(2,5)=40.0, TBCT(3,5)=40.0, 

WBCT(2 ,6)=-0. 197, WBCT(3,6)=-0. 197, 
PBCT(2,6)=3600. , PBCT(3,6)=3600.0, 
FBCT(2,6)=0.0, FBCT(3,6)=0.0, 

IADIX-1, IADIZ-1, 

CV2=63833.0, 


PX(1)=-1.8, PX(2)=0.0, PX(3)=1.8, 

SIZEX(2)=0. 25, 

NXCELL( 1 )=5, NXCELL(2)=5, NXCELT=10, 

PY(1)=0.0, PY(2)=0. 10, NYCELT=1 , 

PZ(1)=1.60, PZ(2)=2.5, NZCELT=10, 

SEND 

Sobs 

N0BS=2 , 

CC(1)=6. 25, CX2(1)=-1 .0, CZ2(1)=-1.0, 

I0H(1)=1 , 


RAL(1)=0. 5, 

CC(2)=6.0, CX2(2)=1.0, 

OPOR(2)=0. 7 , I0H(2)=1 , 
T0BS(1)=0.0, TOBS(2)=100. 
SEND 
$FL 


CZ2(2)=1.0, 

0 , 


CZ(2)=-5.0, 


PRESI=2500.0, 

SEND 

$BF 

SEND 

STEMP 

NTMP=1, ITDIS(1 )=1 , 

TCC(l)=-338.0, TCZ=235.0, TZH(1)=2.5, 
SEND ’ 

SMOTN 
SEND 
SGRAFIC 

NVPLTS=1 , JV1(1)=2, JV2(1)=2, 

NCPLTS=1 , K0NTYP(1)=5, 

JC1(1)=2, JC2(1)=2, 

SEND 

SPARTS 

SEND 


TZL(1)=1.60, 
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Table 1: Major Input and Output Parameters 
for FLOW-3D Tank Pressurization Runs 


Inputs 

Outputs 

Time Step 

Cell Temperature 

Convergence Criteria 

Cell Pressure 

Gravitational Level 

Cell Gas Velocity 

Initial Ullage Temperature 

Average Ullage Pressure 

Initial Pressure 

Fluid Density 

Heat Transfer to Liquid 
Heat Transfer to Wall 
Inlet Gas Velocity 
Inlet Gas Temperature 
Gas Specific Heat 
Mesh Definition 
Obstacle Definition 

Dynamic Viscosity 
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Diffuser 



5 ft Diameter Spherical Tank 

FIGURE 1. - FL0W-3D MODEL OF K-SITE TANK PRESSUR- 
IZATION. 


Low Contour -109.3 R 

High Contour-376. 4 R G — 32.2 ft/sec 2 



Low Contour= 107.9 R 

High Contour =376.1 R G = 3.22 ft/sec' 



Low Contour= 107.4 R 

High Oonlour=376-4 R G = 0.322 ft/sec 2 


Low Contour= 103.3 R 

High Contour =374. 6 R G = 0.0 ft/sec 2 




FIGURE 2, - EFFECT OF GRAVITY ON TEMPERATURE CONTOURS. 55 PERCENT ULLAGE. P = 17 4 psi e 
dt = 0.005 SEC, 2A SEC PRESSURIZATION. 1 






Low Contour - 10ft. 7 It 

High Contour = 3 70.4 [I G 0.322 fl/sec 1 * 


Low Contour 911.30 It 

High Contour- 372.0 It G = 0.0 ft/sec 2 




FIGURE 3. - EFFECT OF GRAVITY ON TEMPERATURE CONTOURS, 28 PERCEN1 ULLAGE, Pj = 17. A psi, e - .007, 
dl = .005 sec, 24 SEC PRESSURIZATION. 


Low Contour= 149.1 R 
High Contour = 626.6 R 


G = 32.2 ft/sec 2 



Low Contour=109.3 R P — 17 4 Dsia. 

High Contour=376.4 R 




Low Contours 183.4 R ^ 

High Contour=503.5R G = 0.0 ft/sec“ 



FIGURE M. - EFFECT OF GRAVITY ON TEMPERATURE CON- 
TOURS, 8.5 PERCENT ULLAGE, Pj = 17.4 PSI, e = .05, 
dt = .003 SEC, 15 SEC PRESSURIZATION. 


Low Contour=313.0 R P, = 1 .] psia 

High Contour=338.5 R r 



FIGURE 5. - EFFECT OF INITIAL PRESSURE ON 
TEMPERATURE CONTOURS, 55 PERCENT ULLAGE, 

G » 32.2 FT/SEC 2 , e = .02, dt = .005 SEC, 
24 SEC PRESSURIZATION. 
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FIGURE 6. - EFFECT OF INITIAL PRESSURE ON 
VELOCITY PROFILE, 55 PERCENT ULLAGE, G = 
32.2 ft/sec 2 , e = .02, dt = .005 sec, 

24 SEC PRESSURIZATION. 



FIGURE 8. - EFFECT OF INITIAL PRESSURE ON VELOC- 
ITY PROFILE, 28 PERCENT ULLAGE, Pj = 17.4 psi, 
e = .007, dt = .005 sec, 24 sec PRESSURIZATION. 


Low Contoux= 108.8 R 

High Cootour=373.4 R P, ~ 17.4 psia 



Low Contour=305.7 R 

High Contour=330.4 R P, = 1.1 psia 



FIGURE 7. - EFFECT OF INITIAL PRESSURE ON TEM- 
PERATURE CONTOURS, 28 PERCENT ULLAGE, G * 
32.2 ft/sec 2 , e = .007, dt = .005 sec, 24 sec 
PRESSURIZATION. 



Low Contour=491,4 R 

High Contour=534.3 ft P, = 1.1 psia 



FIGURE 9. - EFFECT OF INITIAL PRESSURE ON TEMPERA- 
TURE CONTOURS, 8.5 PERCENT ULLAGE, G = 32.2 FT/ 
sec 2 , e = .05. dt = .003 sec, 15 sec PRESSURIZA- 
TION. 
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ORIGINAL PAGE IS 

OF POOR QUALITY 





Pressure, psi 


. t , T, = 250 R at top of tank 

Low Contour = 109.3 R K irk 

High Contour=376.4 r T. = -10 R at liquid surface 


Low Contours 78. 09 R 
High Contour=340.8 R 



FIGURE 10. - EFFECT OF INITIAL TEMPERATURE 
ON TEMPERATURE CONTOURS, 55 PERCENT ULL- 
AGE, Pj * 17.4 psi, G = 32.2 ft/sec 2 , 
e = .02, 24 SEC PRESSURIZATION. 



Time, sec 

FIGURE 11. - RAMP PROFILES USED FOR TEMPERATURE PROFILE 
COMPARISON, 55 PERCENT ULLAGE, Pj = 17 . 4 PS l , 24 sec 
PRESSURIZATION. 


Low Coruour= 109.3 R 
High Coatour=376.4 R 


Ramp Profile 1 



H 


Low Contour^ 113.7 R 
High Contour=393.9 R 


Ramp Profile 2 



FIGURE 12. - EFFECT OF RAMP RATE ON TEMPERA- 
TURE CONTOURS, 55 PERCENT ULLAGE, Pj = 17.4 
PSI, G = 32.2 FT/SEC 2 , e = .02. dt - .005 
SEC. 24 SEC PRESSURIZATION. 
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ORIGINAL PAGE IS 
Of POOR QUALITY 




Q wall — 40 ft-Ili/soc, 



l.uw Contour- 102.4 >| 
High Contours II 71.0 It 


Qli((=4() ft-llj/sui: 




Low Contour i-y 5. 6 5 II Qwail=8U ft-lb/sec, Qliu=<M ft-lb/sec. 

High Contour =369. 5 ft 



FIGURE 13. - tFFECT OP HtAT TRANSFER ON TEMPERATURE CONTOURS, 55 PERCEN1 II Lt AGE , P { - 1/.4 psj, 
G - 32.2 H/sec 2 , e - .02, dt = .005 sec, 24 sec PRESSURIZAIION. 




Low Contour = 108.4 R Qwall=40 ft-lb/sec, Qliq=60 ft-lb/sec 
High Contour =553. 5 R 



FIGURE 14. - EFFECT OF HEAT TRANSFER ON TEMPERATURE 
CONTOURS, 8.5 PERCENT ULLAGE, P, = 17.4 psi, G = 
32,2 ft/sec 2 , e = .05, dt = .003 sec, 15 sec PRES- 
SURIZATION. 


□ 


• □ 


15 sec Pressurization 


O Adiabatic Sokrtton 
□ Qwalt*40 It-lb/'B, Qliq=60 ft-lb/a 
• SLURP Solution 


• o 


& ,° 


o 


• £ 

* o 


• □ 


Distance from Top of Tank, ft 

FIGURE 15. - COMPARISON OF CENTERLINE TEMPERATURES 
USING VARIOUS HEAT TRANSFER RATES, 8.5 PERCENT ULL- 
AGE, e = .07, dt = .005 SEC, Pj = 17.4 PSI. 
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Temperature, 


Low Contours 109.3 R 
High Contour- 376.4 R 


e=.02, dt=.005 sec 


Low Contours 11 7.1 R 
High Contour=376.2 R 


e=.002, dt=.0005 



FIGURE 16. - EFFECT OF CONVERGENCE CRITERIA ON TEMPERATURE CONTOURS, 55 PERCENT ULLAGE, 
Pj = 17.4 psi, G - 32.2 ft/sec 2 , 24 SEC PRESSURIZATION. 



Distance from Top of Tank, ft 

FIGURE 17. - COMPARISON OF CENTERLINE TEMPERATURES 
USING VARIOUS CONVERGENCE CRITERIA, 55 PERCENT 
ULLAGE, G = 32.2 FT/sec 2 , Pj = 17.4 PSI, 24 SEC 
PRESSURIZATION. 










Low Contour 109.3 R ioo cells (10 x 10 Mesh) 

High Contour=376.4 R J 



Low Contour=87,68 R 400 cells (20 x 20 Mesh) 
High Contour=375.8 R 



FIGURE 20. - EFFECT OF MESH SIZE ON TEMPERA- 
TURE CONTOURS, 55 PERCENT ULLAGE, Pj = 

17.4 sec, G = 32.2 ft/sec 2 . e = .02, dt - 
.005 SEC, 24 SEC PRESSURIZATION. 


Low Contour- 103.3 R 100 cells (10 x 10 Mesh) 
High Contour=374 6 R 



FIGURE 21. - EFFECT OF MESH SIZE ON TEMPERA- 
TURE CONTOURS, 55 PERCENT ULLAGE, Pj - 17.4 
psi, G = 0.0 FT/SEC 2 , e = .02, dt = ,005 
SEC. 24 SEC PRESSURIZATION. 
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